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SUMMARY 

Analytical methods for the determination of local values of 
outside and inside heat- transfer coefficients and effective gas tem- 
peratures from turbine-blade- temperature measurements were developed. 
Temperature-distribution equations are derived for typical turbine- 
blade configurations at the central section when the wall thickness 
is uniform, at the leading-edge section, and at the trailing-edge 
section. The equations have the same general form for all sections 
of the blade, but the blade configuration at the location where the 
data measurements are taken affects the evaluation of some terms in 
the general equation. 

Procedures for applying these analytical methods to experimen- 
tally measured blade-metal temperatures* are presented. Data are 
presented for the leading and trailing edges of a symmetrical water- 
cooled blade to illustrate the validity of the methods for those 
portions of the blade. 

In addition to the application to turbine blades, the methods 
can be applied to any heat-transfer apparatus having a profile that 
can be approximated by the shapes discussed. 


INTRODUCTION 

The development of a sound basis for the design of cooled tur- 
bines is dependent on knowledge of the inside and outside heat- 
transfer coefficients for turbine blades. Average heat- transfer 
coefficients can be determined for any particular blade configuration, 
but it appears infeasible at present to obtain a correlation of average 
coefficients that is suitable for all blade configurations at all tem- 
perature ratios and all velocity distributions. The fundamental 
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boundary- layer and heat- transfer theories are therefore being studied 
at the NACA Lewis laboratory as a means of predicting local heat- 
transfer coefficients for any given set of flow conditions. These 
studies will eliminate the necessity for experimental work on every 
conceivable blade shape. 

In order to verify and to extend these theoretical studies, 
local heat-transfer coefficients around the periphery of a limited 
number of turbine blades must be determined over a wide range of 
flow conditions. 

Three methods were proposed for obtaining these coefficients: 

(1) boundary- layer surveys using a hot-wire anemometer, (2) boundary- 
layer surveys using an optical interferometer, and (3) analytical 
solutions using measured turbine-blade temperatures. Considerable 
progress has been made in the development of operating techniques 
for the hot-wire anemometer and the interferometer, but some prob- 
lems must still be solved. The third method is probably the best 
means presently available. A development of analytical solutions 
is presented herein for calculating local values of outside- and 
inside- convection heat- transfer coefficients from experimental data 
obtained at steady-state conditions from the measured temperatures 
of turbine blades having a known thermal conductivity. 


METHODS OF ANALYZING HEAT-TRANSFER DATA TO OBTAIN LOCAL 
CONVECTION COEFFICIENTS AND EFFECTIVE GAS TEMPERATURES 

The quantity of heat transferred per unit area by convection 
from a gas stream to a solid surface can be expressed as the product 
of the surface heat-transfer coefficient and the effective tempera- 
ture difference between the gas and the surface. For one- dimensional 
heat flow, this same quantity of heat will be conducted through the 
wall of the solid and can be evaluated by multiplying the thermal 
conductivity of the wall material by the temperature gradient in the 
wall. If the inside of the wall is convection-cooled by a fluid, the 
heat transferred can also be expressed as the product of the inside- 
surface heat-transfer coefficient and the effective temperature dif- 
ference between the inside surface and the coolant. By using this 
reasoning, a heat balance can be set up so that the inside and out- 
side heat-transfer coefficients can be calculated from a known tem- 
perature difference between two known locations in the wall having 
a known thermal conductivity. 
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The temperature-distribution equations for any vail configura- 
tion where the heat flow -is one dimensional can be reduced to the 
following form (Symbols are defined in appendix A and equations are 
derived in appendix B.): 

t x - ty = r(t e - t x ) (l) 

and 



where 

t x and ty. wall temperatures at two specific locations lying in 
direction of heat flow 

tg wall temperature any place in line of heat flow 

(usually taken as either t x or ty) 

T proportionality factor that is function of wall con- 

figuration, location of temperature measurements 
and ty, and ratio of outside heat-transfer 
coefficient h^. to wall thermal conductivity kg; 

can be mathematically expressed and experimentally 
determined 

* / 

£ and q proportionality factors that are functions of wall con- 
figuration, location of temperature measurement tg, 
wall thermal conductivity, and ratio of outside heat- 
transfer coefficient to wall thermal conductivity; 
can be mathematically expressed and evaluated by use 
of ratio ht/kg calculated from experimental deter- 
mination of r 

The mathematical expressions for T, £, and q are relatively 
simple for a simple shape such as a wall of uniform thickness, but 
they become more complicated for walls defined by concentric circles 
or for rectangular- or trapezoidal-shaped fins. Mathematical expres- 
sions for shapes other than those mentioned have not been derived. 

For a given configuration and thermocouple location, T, k;g£, and 
kgq can be mathematically evaluated as functions of h-t/kg and 
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plotted over the range of ht/kg that can be expected in experi- 
mental investigations. This fact greatly simplifies the use of 
these quantities. 

From equation (1) it can be seen that if t x , ty, and t e 
are known from experiment, the value of T is easily calculated. 

The quantity h^/kg can then be obtained for a mathematically deter- 
mined curve of T as a function of h t /kg.- For a known value of the 
thermal conductivity kg, the outside heat-transfer coefficient h t 
is then also known. The values of ht/k;g and kg can then be used 
to evaluate £ and tj, and if the coolant temperature t c is known 

from experiment, the inside heat-transfer coefficient h^ can be 
calculated from equation (2). 

Frequently the effective gas temperature t e is unknown and 
will also have to be evaluated from experimental investigations. 

In order to evaluate t e , a series of experimental data points must 
be obtained for a constant outside heat- transfer coefficient h^. 

From equation (1), if t x - ty is plotted against t x for a constant 
value of T, the intercept on the t x -axis will be t e (this is the 
case where t x - t y = 0) and the slope of the plotted line will be 
-r. By using this method it is possible to determine both h-t (from 
— r) and tg . 

It is often more convenient to write equation (l) in the form 



where tg is an observed gas temperature, and 

n = ^ (4) 

t g 

From equation (3), if (t x "ty)/tg is plotted against t x /tg, 
the intercept on the t x /tg-axis will be Q and the slope -will 
again be -T. The effective gas temperature is then calculated from 
the value of Q by use of equation (4), and the outside heat- transfer 
coefficient is determined from the slope -T. 

Sometimes measuring the temperature at two points in the wall 
is impossible. If only one wall temperature tg is known, the 
temperature-distribution equation for determining the outside heat- 
transfer coefficient can be written 
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- t c = r(te " tc) (5) 



For the cases vhere the temperature gradient in the metal could 
be measured by means of two thermocouples, the term T was a function 
of ht/kB and the blade dimensions. For the case with one wall- 
temperature measurement, however, T is also a function of the inside 
heat-transfer coefficient h^. Consequently, the method of applying 
equation (5) or (6) to experimental data is different from the method 
used for equations (l) and (3). The procedure is rather complex and 
has been utilized for a uniform wall thickness only; it is explained 
in detail in the section entitled "Application of Heat-Transfer Equa- 
tions to Experimental Data." 

The temperature-distribution equations used for all these ana- 
lytical methods are based on a total outside- surface heat-transfer 
coefficient that is a combination radiation and convection coefficient. 
Corrections can be made to this total heat-transfer coefficient for 
radiant-heat transfer to obtain a true convection coefficient. 

These methods of calculating heat-transfer coefficients and 
effective gas temperatures can be applied to any apparatus where the 
heat flow is one dimensional and the heat transfer is through a shape 
that can be approximated by the simple shapes discussed herein. The 
central section (with uniform wall thickness), the leading-edge sec- 
tion, and the trailing-edge section of most turbine blades can be 
approximated by simple shapes so that these methods of analysis can 
be used to determine local values of inside and outside heat- transfer 
coefficients and effective gas temperatures. 


Central Section of Turbine Blade 

Equations are presented for blades having a uniform wall thick- 
ness between the leading-edge and trailing-edge sections for plain 
hollow blades as shown in figures 1(a) and 1(b). Analyses for the 
central sections of other blade configurations have not been verified 
and therefore will not be presented here. The derivations for all 
equations are given in appendix B. 

Case for measurable wall-temperature gradient. - In order to 
use this method, the wall temperature must be known at two points 
located in the line of heat flow. 
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The heat flow through the blade wall is assumed normal to the 
blade surface, and the surface curvature is assumed small enough 
that the wall can be treated as a flat plate. With the substitution 
of temperatures tq and t2 at the locations shown by small circles 
on figure 1(a), equations (l) and (3) take the form 


and 


where 


*1 - t 2 = T(t e - tq) 



r = 




(6i - & 2 ) 


1 + ^ (6 - 6i) 

k B 


The dimensions bq, &2, and 
The values of £ and q 


6 are shown in figure 1(a). 
in equation (2) are 


and 



(7) 

( 8 ) 


(9) 


( 10 ) 

(ID 


Case with one temperature measurement in blade wall. - The 
temperature-distribution equation for the case with one thermocouple 
in a wall of uniform thickness can be written 


or 


where 


^B “ t c ~ ^(^e 


(5) 

t B - t c _ r (b 

t rp V 

. 

t a ) 

(6) 


r = hj-Ckg + &h ± ) 

k B (h t + hi) + 6 h t hi 

The dimensions 5 and b are shown in figure l(b). 


( 12 ) 
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Leading-Edge Section of Turbine Blade 

The heat is assumed to be transmitted from the gas stream 
through the blade at the leading edge to the coolant along a sector 
having an included angle d0, as shown in the cross-sectional view 
of the leading-edge section in figure 1(c), where thermocouples 3 
and 4, indicated by small circles, lie in a direct line between the 
stagnation point on the blade surface and the coolant passage. The 
equation for the temperature distribution at the leading edge can 
be written 
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t3 - t4 = r(t e - t3 ) 


(18) 


or 


h : h , r/n - 


(19) 


and 


r = 


r o h t . r 3 

~T~ loge 77 

k B r 4 


1 + 


r o h t r o 


r 3 


k B 


l°g e 


where 

r radii as shown with subscripts on figure 1(c) 
The values of £ and tj in equation (2) are 


i -cr 


r i ( k B 


k B W~o + 1086 


and 


( 20 ) 


( 21 ) 


_ r i / k B 


B = ~ 


k B \ h t r o 


+ log e -2 


( 22 ) 


Trailing-Edge Section of Turbine Blade 

A cross section of the trailing-edge section of most turbine 
blades can be very closely approximated by a trapezoid, a rectangle, 
or a combination of trapezoids and rectangles. Temperature distribu- 
tions have been determined for rectangular and trapezoidal cross 
sections where the heat flow is assumed to be one dimensional. Der- 
ivations are given in appendix B for trailing-edge sections composed 
of the following shapes: 

1. Trapezoidal 

2. Bectangular 
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3. Combination of one trapezoid and one rectangle 

4. Combination of tvo trapezoids 

Equations for T, £, and tj are given for each of the shapes. 


Trapezoidal trailing-edge section. - By referring to the trap- 
ezoidal trailing-edge section on figure l(d) for thermocouple loca- 
tions and dimensions, equations (l) and (3) can be written 

ts ~ = r*(te - t5 ) (23) 


and 


where 




(24) 


(25) 


The values of £ and T) in equation (2) are 


l - f (26) 

and 

B = f (27) 

where 

N = |~Hi(i|i) Jo(i|) + iJl(Ul) iHo(i|)] ( 28 ) 

I = iJ-^(it^) H^i^) - H^(i£^) iJ^(i|2) (29) 


iz 

2B 2 k B 


|”ffl(i!l) Jo(i|2) + iJ l( i ll) iH o( i l2^| 


E 


(30) 
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| = 2B 


/ It _ 

\ T ( 

1 - tan a\ 

\|V y + J T ’ 

l) + T 1 [ 

2 tan a ) 


B 2 = 


ks sin a 


and 


a = tan 


-1 T 3 " T 1 

2L 


n 6 

£2 

7, 


V V 


evaluated for y = y 3 

evaluated for 7 = 7q 

tcT-| 

evaluated for y + — ± = 0 

4 

evaluated for y = L 

shown on trailing-edge section sketch for trap- 
ezoidal shape (fig. 1(d)) 


Where there are no subscripts on N and they can be 

evaluated for either thermocouple in the trailing-edge section. The 
resulting value of h^ in equation (2) should be the same in either 
case. 


Rectangular trailing-edge section. - For the rectangular trailing 
edge section as illustrated by the trailing-edge section shown in com- 
bination with a trapezoidal section in figure 1(e), equations (l) and 
(3) are written 

^5 " ^7 = “ ^ 5 ) ( 3 l ) 

and 

ts ' , r(b . (32) 

and 

r = cosh H 7 ? + f T i) - cosh <p( 7 5 + f T i) ' (33) 

COSh <p(y 5 + * tJ 
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The values of £ and q in equation (2) are 

cosh qp(y + f T i) 


t = 


k B cp sinh cp^L 2 + J T 1 ) 


(34) 


and 


cosh <p(l 2 h- l 
k B Cp sinh op(l 2 + f T i) 


(35) 


where 


<p = 


|2ht 

*B T 1 


The dimensions y, T, and L are shown in figure 1(e). 

Combination trapezoidal and rectangular trailing-edge section. 
For rectangular portion d of the trailing-edge section shown on 
figure l(e), the value of T is the same as in equation (33). 


r = r d = 


_ cosh cp d ( y ? + | Tl ) - cosh cp d (^y 5 + f 


cosh CPd(y5 + J T l) 


(33a) 


For the trapezoidal portion b, 

t 7 " t 6 = ^b(t e - t 7 ) 


(36) 


or 


1 7 “ t 6 


■ r A n - i 


(37) 


and 


r b ' 


J 0 (l 5b,6> 


J 0^b,7^ + dK lH 0^b,6^ " lH 0^b,7^ 


+ f§ i ^o( i ^b,7) 


( 38 ) 
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where 


2B b 21J l 

K,1 

) + *b,l *0 

L [tanh cp d 

( L 2 + 4 T l) 

/o' 

^,i) 

2B b 2 H l( 

,^b,l) 

- *b,i *d[ 

tanh cp^ 

L 2 + f T l)_ 

iE o( i ^b,l) 


(39) 


h ’ 2 h\jy * T i(: 


1-tan a-^ 
2 tan 


-1 t 3 ‘ T 1 
% = tan ± _ ± 

2L X 

and is evaluated for y = 0 . ' The subscript b refers to 

portion b of trailing- edge section. 

The value of in equation (38) is a function of both 

ht,b/kB and kt,d/ k B* For most accurate solution to the equa- 
tion, it is necessary to evaluate dz/dEC in equation (39) using the 
value of htjd/kB obtained for the rectangular portion of the trail- 
ing edge by use of equation (33a) . It is doubtful that this much 
trouble is warranted, however, because the terms in equation (38) 
involving dz/dK are small compared with the other terms in the 
equation so that the error caused by assuming ht,d = ht,b would 
be negligible. By using this assumption, 1^ can be evaluated as 
a function of h-t^ only. 

The method for determining the value of h^ at the trailing- 
edge where the section is composed of two different portions is so 
insensitive to the metal-temperature measurements that the analytical 
solution is not believed to be worthwhile and therefore is not 
presented. The determination of hi is much more accurate at the 
central portion of the blade where the shape factor is simpler. The 
evaluation of the outside coefficient ht is considered to be quite 
accurate at the trailing-edge section, however. Because such a small 
portion of the coolant passage is next to the trailing-edge section, 
the evaluation of the inside heat-transfer coefficient at that por- 
tion of the blade is of minor importance. 

Combination of two trapezoidal portions in trailing-edge sec- 
tion. - For trapezoidal portion d* on figure 1(f), 
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i 


or 


and 


where 


t 5 " t 7 = r d' 



(40) 


(41) 


(42) 


N = 


Ed' ,2 

2B d ' 2 k B 


[aifiSa'.i) 


^O^d' ) + iJ l^ i ^d',l^ 





k B sin a^-i 


$d' 



and 


a 


d* 



For trapezoidal portion h' on figure 1(f), 
t7 - t6 = r h’(t e - t7) 
or 



(43) 


( 44 ) 
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J 0^b',6^ " J 0^ i ^b',7^ + “ iB o(^b' ,7)] 


where 


+ Sr ^^b',?) 


dZ' GlJ i( 1 ^b , ,l) + SJ 0^b',l^ 
= GHili5 b .,i) - SiHod^,^) 


n -c ft \ 

G = Jo^dSz) + H^tiidS!) 


B d' Sb',1 ^l^d'.l) ^^^*, 2 ) , , ,., \ 

— h M„.) - 1J i (i| as 2 ) 


'l'-^d',!' 


Also 


£ b’ “ 2B b' ^ / y + T : 


(1 - tan 
2 \ 2 tan a- , 


t = tan 


T - T 
-1 3 '2 


where 



^b',6 

evaluated for 

y = ^6 

tb' ,1 

evaluated for 

y = y 7 

£b' ,1 

evaluated for 

0 

11 

>» 
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2 evaluated for y = L 2 


JtT -i 

1 evaluated for y + — - = 0 

A method of evaluating h.^ for this section of the blade is not 
presented for reasons previously discussed. 


Radiation Corrections 


The effects of radiation have been neglected up to this point 
in the analysis so that the local outside heat-transfer coefficients 
obtained are combination convection and radiation coefficients. In 
most cases, the heat transferred by radiation is relatively small 
(about 3 to 10 percent of the heat transferred by convection for 
uncooled surface temperatures up to 1500° F), but it is still of 
sufficient magnitude to require evaluation. 

The combination convection and radiation coefficient has been 
defined as h t ; if the convection coefficient is designated h G 
and the radiation coefficient, h r , the following equations may be 
written: 


h t = h o + h r 


and 


h-t- = 


Qt 

A B,o^e " V 


Q + Q r 

A B,o^e " V 


(47) 


(48) 


where 


Q-u total heat-flow rate to blade 
Q heat-flow rate to blade by convection 

Q r heat-flow rate to blade by radiation 

The heat transferred by radiation is given in reference 1 as 
Q r = 0.173 A BjC 
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therefore 



The actual value of the convection heat-transfer coefficient h D 
can now he calculated, for each gas temperature from equations (47) 
and (50). 

A comprehensive discussion of radiation is unwarranted here. 

For a precise evaluation of radiation, it is necessary to accurately 
determine the geometry factor F and the emmissivities. A method 
of determining geometry factors by the use of a mechanical integrator 
on large-scale models is suggested by Hottel in reference 2. For 
most applications with cooled turbine blades, the metal surfaces will 
be tarnished and soot covered so that the eramissivity will be high, 
probably ranging from 0.80 to 0.95 for nearly all materials. 


APPLICATION OF HEAT-TBANSFER EQUATIONS TO EXPERIMENTAL DATA 

There are two methods of applying the heat-transfer equations 
to experimental heat-transfer data. The most direct method, which 
is probably also the better method, depends upon a knowledge of the 
effective gas temperature. The, usual method of determining the 
effective gas temperature is calculation by use of a known blade 
recovery factor. A discussion of recovery factors is contained in 
reference 3. When the effective gas temperature, the measured blade 
temperature at two positions in the direction of heat flow, and the 
metal thermal conductivity are known, the outside and inside heat- 
transfer coefficients can be calculated directly. The second method 
can be used if the effective gas temperature is unknown, but data are 
necessary from heat- transfer runs where the heat- transfer coefficient 
on the outside surface of the blade is maintained constant while the 
amount of heat transfer to the blade is varied by allowing the coolant 
flow and either the coolant temperature or the gas temperature to 
change . 
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The second method cannot be used to determine the effective 
gas temperature with as great an accuracy as is possible using a 
local blade recovery factor. The method is therefore best applied 
when other means of determining effective gas temperature fail, such 
as the case where temperature gradients in the gas stream make meas- 
uring the stream temperature at the blade impossible. 

/ 

Case with Known Effective Gas Temperature and 
Measurable Wall-Temperature Gradient 

The effective gas temperature, or the gas temperature effecting 
heat transfer, is defined as the adiabatic surface temperature in 
reference 3; that is, the effective gas temperature is the tempera- 
ture that the surface would assume if it were thermally insulated 
so that there would be no heat transfer. The local effective gas 
temperature can be calculated from the total gas temperature, the 
total pressure, and the local static pressure if the local blade 
recovery factor is known. For simple shapes such as flat plates, 
tubes, and wedges, analytical solutions for the recovery factor are 
available (reference 3). For more complex shapes, such as turbine 
blades, recovery factors can be determined from adiabatic tests of 
the blades. Experimental work at the Lewis laboratory has shown 
that turbine-blade recovery factors are usually quite close to 0.90 
for all Reynolds numbers and for local Mach numbers from 0.4 to 1.0. 
The effective gas temperature is then 

t e = t + A b (T -*t) (52) 

or 



r 

7-1 

1 

t e = T ' 

1 - (1 - A b ) 

7 


1 'I) 



V. 

— N-*- / — 

- 


In actual practice for subsonic flow, the value of the effective 
gas temperature t e will be 

0.98 T < t e < T 

so that in many cases the error in assuming t e = T will be neg- 
ligible. The size of the error will be dependent on the temperature 
difference t e - t x . A more complete discussion of this error will 
be found in the section entitled "Accuracy Considerations." 


\ 
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Equation (l) can be written 


r = IcJa d») 

t e ” z x 

so that T can be evaluated directly from the effective gas tem- 
perature and the measured blade temperatures. For a given config- 
uration and given thermocouple locations, T can be mathematically 
represented as a function of h^/kg so that the heat-transfer coef- 
ficient ht can be easily calculated. The inside heat- transfer 
coefficient h^ can then be evaluated from 



where £ and q are obtained from mathematically determined curves 
of kg£ and k^Ti plotted against ht/kg for the given configura- 
tion and thermocouple location. 


Case with Unknown Effective Gas Temperature and 
Measurable Wall-Temperature Gradient 

The method of application for this case is considerably more 
complex than for the case where the effective gas temperature is 
known. As previously stated, a constant heat- transfer coefficient 
must be maintained on the outside surface of the blade for a series 
of experimental heat-transfer runs in order to utilize this method 
of analysis. If the coolant flow and the coolant temperature are 
allowed to vary. for a series of runs, the outside heat-transfer 
coefficient can be held constant by maintaining the gas flow and 
the gas temperature constant. The values of inside and outside heat- 
transfer coefficients and the effective gas temperature can then be 
obtained rather easily by use of equation (l) 

t x " V = ^e “ (- 1 -) 

which can be solved by use of suitable plots. If t x - ty (ordinate) 
is plotted against t x at a constant value of T (constant ht/kg) 
over a range of blade-metal temperatures, the slope of the resulting 
line will be -T and. the intercept on the t x -axis will be the 
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effective gas temperature t e . The value of the local outside heat- 
transfer coefficient h-£ can now be calculated from the equation 
for r for the section of the blade in question. 

In order to obtain a sufficient number and a suitable range of 
temperatures t x and ty to determine the line in the plot of 
t x - ty against t x , it is necessary to vary both the coolant flow 
and the coolant temperature over as wide a range as possible. 

After determining ht and t e , the value of the local inside 
heat-transfer coefficient hj_ can be determined from equation (2). 

The values of £ and q are evaluated as functions of h^/k-g for 
the section of the blade. 

If no provision is made for varying the coolant temperature, 
the gas temperature must be varied. The outside- surface heat- transfer 
coefficient is apparently affected by the ratio of gas temperature to 
blade temperature; therefore it may be necessary that the weight rate 
of gas flow be varied at the same time the gas temperature is varied 
in order to maintain a constant value of outside-surface heat- transfer 
coefficient. By making heat- transfer runs over a large range of gas 
flow and gas temperature, the effect of the ratio of gas to blade 
temperature can be determined so that data for a constant value of 
the outside-surface heat-transfer coefficient can be obtained for 
the required plots. 

This effect can be evaluated if the variation in average con- 
vection coefficients with temperature ratio is assumed to be the same 
as the variation in local convection coefficients. This assumption 
will be valid if the point of transition from laminar to turbulent 
boundary- layer flow on the outside surface of the blade remains at 
the same place for all temperature ratios. Evidence exists indi- 
cating that the amount of laminarity is a function of temperature 
ratio unless the pressure gradient is strong enough to maintain a 
completely laminar boundary layer. Where the pressure gradient would 
indicate that the boundary layer is either almost completely laminar 
or almost completely turbulent, the following procedure can be used 
for evaluating the effect of temperature ratio on the outside heat- 
transfer coefficient: 


1. The average outside heat- transfer coefficient h^. av can 
> be calculated from the equation 


c p,c w c At c 


B,av 


y 


f 
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where 

t e = 0.98 tg . (assumed value of ft = 0.98 as an approximation) 

Then ht,av can he plotted against the weight rate of gas flow Wg 
for various measured gas temperatures tg as shown in figure 2(a), 
Plot I. 

2. From Plot I, a cross plot can he made as shown in fig- 
ure 2(h), Plot II, to obtain tg as a function of Wg for various 
values of h-t^g-y. 

3. Local hlade temperature tg and coolant temperature t c can 

be plotted against Vg for the various measured gas temperatures at 
a constant coolant flow, as shown in figures 2(c) and 2(d), Plots III 
and IV, respectively. 

From Plots II, III, and IV, it is possible to correlate the 
weight rate of gas flow and the temperatures of the hlade, coolant, 
and gas for given constant values of h-^ aV and coolant flow. These 
plots are required for obtaining temperatures to use in additional 
plots to determine outside and inside heat-transfer coefficients and 
effective gas temperatures. 

By use of equation (3) 



it can he seen that if (t x - ty)/tg (ordinate) is plotted against 
t x /tg for a wide range of gas temperatures and coolant flows at a 
constant value of ht,av> 'the slope of the resulting line will he 
-r and the intercept on the t x /tg-axis will he ft. The values 
of ty, t x , and tg are obtained from graphs like Plots II and 
III of figure 2. 

The value of the local outside heat-transfer coefficient h^ 
can be obtained as before from the equation for T for the section 
of the blade in question, and the effective gas temperature is cal- 
culated from 


t e = ftt 


(4a) 
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In order to simplify the evaluation of h-t from a known value 
of T, it is convenient to plot T as a function of ht/kg for 
the particular section of the blade, as previously explained. 

The value of the local inside heat-transfer coefficient hi 
can be determined from equation (2) where the value of t c is 
obtained from Plot IV of figure 2, ht and t e are determined from 
equation (3), and the values of £ and tj are evaluated as functions 
of ht/kg for the blade section. 


Case with Known Effective Gas Temperature and 
One Temperature Measurement in Blade Wall 
Equation (5) can be written 



(5a) 


so that T can be evaluated directly from the effective gas tem- 
perature, the coolant temperature, and the blade temperature. Because 
T is a function of both hi and ht/kp, a plot is necessary to 
determine the outside heat-transfer coefficient h t . The data for 
this plot must be obtained from a series of experiments where the 
outside heat-transfer coefficient is maintained constant while the 
coolant flow is varied, which varies the inside heat-transfer coef- 
ficient. The value of T is then calculated from equation (5a) for 
each experimental point. The mathematical equation for T is 


p ht(kp + 8 hi) 

k B (ht + h t ) + 6 h.^ 

which can be written 

1 - r _ r /_c_ + SC\ _ 5C 

w c n \ h t k B/" k B 

where hi was replaced by Cw c n . The exponent n can be evaluated 
from calculations of the product of the average inside coefficient 
and the coolant-passage area for the variable coolant-flow runs. 

(The area need not be evaluated as it is constant.) 


( 12 ) 

(15) 
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A,h, 


'P.e 


w_At, 


1 1 ,av 


(53) 


J B,av 


If A^h^ av is plotted against w c on logarithmic coordinates, 
the slope of the resulting line is the exponent n. 


It can he seen from equation (15) that if (1 - r)/w c n is plotted 
against V for variable coolant-flov runs, the intercept on the 
T-axis of the resulting straight line is T' (for (l - r)/w c n = 0). 
Substitution of this value of T' into equation (15) yields 



(54) 


from which the value of h^ can be calculated. 

The local inside heat-transfer coefficient hj_ is then cal- 
culated from equation (2) using equations (16) and (17) to evaluate 
J and T) , respectively. 

This method of determining heat-transfer coefficients using only 
one blade-metal- temperature measurement, however, is not completely 
satisfactory in all cases. The method works best for blades made of 
metals having low thermal conductivities. When the blade metal is 
thin and the metal thermal conductivity is high, obtaining an accu- 
rate evaluation of the outside coefficient h-t becomes difficult, 
although the ratio of inside to outside coefficients can be quite 
accurately evaluated. If the local inside coefficient hj_ can be 
obtained from measurements in the coolant passage, the value of h-^ 
can then quite easily be obtained from equations (5a) and (12). 


Case with Unknown Effective Gas Temperature and 

One Temperature Measurement in Blade Wall 

When the effective gas temperature is unknown, the experimental 
values of local heat-transfer coefficient and effective gas tempera- 
ture are obtained from a graphical solution of equation (6) 

r ( n - l i) < 6 > 



i 
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as illustrated in figure 3. In Plot I (fig. 3(a)), (tg - t c )/tg 
is plotted against t c /tg for various constant values of coolant 
flov (and thus constant hp), where t c , tg, and tg are obtained 
from plots like those in figure 2 for a constant value of h^ av . 

The slope of each line is -T and the intercept of all lines 
on the t c /tg-axis is Cl. The effective gas temperature t e can 

now be calculated from 


t e - Qtg (4 a ) 

The graphical solution to equation (15) is illustrated in Plot II 
(fig. 3(b)). The values of T for corresponding values of w c are 
obtained from Plot I (fig. 3(a)). The straight line representing the 
plotted values of T against (i -r)/w c n on Plot II can be extended 
to intercept the T — axis at a point denoted by I”, as previously 
explained. The value of ht can then be calculated from equa- 
tion (54). The local inside heat-transfer coefficient hi is then 
calculated from equation (2) in the same manner as previously 
explained. Again, this method is not completely satisfactory for 
all cases, particularly when the metal thermal conductivity is high. 


ACCURACY CONSIDERATIONS 

In the experimental determination of heat- transfer coefficients 
from cooled turbine blades, sources of error exist in the experi- 
mental measurements and in the method of analysis that must be 
minimized. These possible sources of error are: 

1. Assumption of one -dimensional heat transfer 

2. Effect of approximations to blade configuration 

3. Effect of variations of thermal conductivity in direction 

of heat flow 

4. Effect of variation in heat- transfer coefficient along blade 

in chordwise direction 

5. Effect of temperature gradient in trailing-edge section 

6. Accuracy required in locating thermocouples 
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7. Effect of thermocouple hole on heat-flow path 

8. Accuracy of temperature measurement 

Each of the preceding items is discussed in detail. 


Assumption of One-Dimensional Heat Transfer 

Errors in the assumption of one- dimensional flow are negligible 
except near the rim on a cooled turbine, as shown in references 4 and 
5. The lower the metal thermal conductivity, the smaller the errors 
involved. The metals used in nearljr all turbines do have a low con- 
ductivity. In static-cascade work, some means of thermally insulat- 
ing the ends of the blades is usually used so that spanwise tempera- 
ture gradients are minimized even at the blade ends. 

The temperature distribution through the cross section of a 
typical water-cooled turbine blade was computed by relaxation methods 
(reference 6) for two metal thermal conductivities and for both a 
constant and a variable outside heat-transfer coefficient. In all 
cases, the isothermal lines showed that the heat flow was very close 
to being one dimensional at the leading and trailing edges. It was 
also found that the temperature distribution in the trailing-edge 
section of the blade could be accurately determined by a calcula- 
tion of the temperature distribution through a wedge of comparable 
dimens ions . 

At the central section (midchord) of a blade with a uniform 
wall thickness, the heat flow will be one dimensional except at the 
ends of the coolant passage next to the leading- and trailing-edge 
sections. At such locations, where the flow is known to be two or 
three dimensional, these analytical methods should not be used. 

The assumption of one- dimensional heat flow is valid when 
reasonable care is used in locating thermocouples in regions away 
from end effects. 


Effect of Approximations to Blade Configuration 

As previously mentioned, the temperature distribution in the 
trailing-edge section can be accurately determined by a calculation 
of the temperature distribution through a wedge (trapezoidal trailing- 
edge section) of comparable dimensions (reference 6). Conversely, if 



1347 


NACA EM E50F09 


25 


the temperature distribution is known (defined herein by two tem- 
peratures at two known locations), the outside heat- transfer coef- 
ficient can be accurately determined. The same assumption of one- 
dimensional heat flow is used for a rectangular trailing edge, and 
therefore the same accuracy should be obtained. Trailing-edge cross 
sections can usually be approximated very closely by some combination 
of rectangles and trapezoids, so that errors due to shape approxima- 
tions can be made negligible. 

At the leading edge of the blade, the accuracy will depend upon 
the configuration. If the leading edge is well rounded, relaxation 
calculations indicate that the shape approximation of two concentric 
arcs will give results accurate to about 10 percent. The applica- 
tion of the leading-edge analysis is less accurate than the applica- 
tion of analyses for the central and trailing-edge sections because 
the shape approximation is not completely valid. 

The assumption that the central portion of the blade can be 
approximated by a flat plate is accurate for large radii of curva- 
ture. For small radii of curvature at the central portion of the 
blade, more accurate results can be obtained by use of the leading- 
edge equations for arcs of concentric circles. The leading-edge 
equation approaches the central-section equation as the radius 
"* • approaches infinity. The approximate percentage error involved in 

using the flat-plate equation for determining T can be determined 
from the expression 


E rr or ( pe r c ent ) 



(55) 


where 

6 wall thickness 

r Q radius of curvature of outside surface of blade 

For a given value of T, the calculated heat-transfer coefficient h^ 
is always higher for the flat-plate equation. 


Effect of Variations of Thermal Conductivity 
in Direction of Heat Flow 

The analytical methods were derived with the assumption that 
the section of the blade in question was of a material having a uni- 
form thermal conductivity. The thermal conductivity of a material 
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is a function of its temperature. There must be a temperature gra- 
dient to have heat flow; therefore the thermal conductivity varies 
in the direction of heat flow. Analyses have been made to evaluate 
the effect of variable thermal conductivity caused by temperature 
gradients. For a typical high-temperature-metal alloy with the 
thermal conductivity given by 

k = 2.763 X 10" 6 T + 1.246 X 10“ 3 (56) 

the temperature distribution was calculated for temperature differ- 
ences as high as 130° F between two points in the metal. The error 
caused by using an average metal temperature for evaluating the 
thermal conductivity was negligible. 

If the blade is laminated or coated on the inside or outside 
so that a nonuniform thermal conductivity is caused by the use of 
different materials, a great deal of care will be required in inter- 
preting the results of the analysis. The heat flow can still be 
measured as long as the temperature measurements are made in a piece 
of material of constant thermal conductivity, but the heat-transfer 
coefficients obtained from the analysis will contain the effects of 
convection, radiation, and conduction through the material having a 
thermal conductivity different from the material in which the tem- 
perature measurements are made. 


Effect of Variation in Heat-Transfer Coefficient 
along Blade in Chordwise Direction 

At the leading edge of the blade where the heat- transfer coef- 
ficient changes rapidly with the distance from the stagnation point, 
the determination of the outside heat-transfer coefficient will be 
inaccurate unless the thermocouples used to measure the temperature 
gradient lie in a direct line between the stagnation point and the 
coolant passage. If the leading-edge section is long and relatively 
sharp, the results obtained from the analysis will be of very doubtful 
quality. The analysis method is based on a short, well-rounded 
leading-edge section. 

At the central and trailing-edge sections of the blade, the 
variation in heat-transfer coefficient in a chordwise direction is 
usually gradual so that its effect will be quite small. At the cen- 
tral portion of the blade the effect will be negligible and at the 
trailing-edge section the calculated heat-transfer coefficient will 
be an average value for the blade surface between two thermocouple 
stations. The closer these stations are together, the smaller will 
be the effect of variation in heat-transfer coefficient. 
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Effect of Temperature Gradient in Trailing-Edge Section 

A decrease occurs in the magnitude of local heat- transfer coef- 
ficients as the surface temperature of a flat plate increases in the 
direction of fluid flow for a laminar boundary layer (reference 7). 
This decrease in heat- transfer coefficient is caused by the formation 
of a cool boundary layer. that serves to insulate the surface from the 
gas stream and thus decreases the amount of heat transfer. 

This decrease in heat transfer can be defined, by a decrease in 
effective gas temperature rather than a decrease in heat- transfer 
coefficient; but it is believed that defining the effective gas tem- 
perature using a recovery factor is a better procedure. The tem- 
perature gradient should then cause a decrease in heat-transfer 
coefficients along the trailing-edge section. 

By using the analytical method incorporating a known effective 
gas temperature, this effect of decreasing heat-transfer coeffi- 
cients can be experimentally determined by placing the thermocouple 
stations close together. By using the method where the effective 
gas temperature and the heat-transfer coefficients are simultane- 
ously determined, the effective gas temperature would probably be 
abnormally low because of the effect of the cooled boundary layer, 
and the heat-transfer coefficient would be higher than for the case 
using a known effective gas temperature. The calculated rate of 
heat transfer to the blade would be the same in either case, but 
care should be exercised in evaluating the results obtained by the 
two different methods. 


Accuracy Required in Locating Thermocouples 

The distance between any two thermocouples used for measuring 
a temperature gradient should be known as accurately as possible 
because the errors in the resulting data are in direct proportion 
to the error in the measurement of the distance between the loca- 
tions where the temperature measurements are taken. Because ther- 
mocouples are usually located inside drilled holes, it is advanta- 
geous to have relatively shallow holes in order to reduce the amount 
of drift in the drilling operation as much as possible. 

The location of the thermocouple junction in the hole is also 
of prime importance. The thermocouple junction should therefore be 
as small as possible and the junction itself should be in a hole of 
t reduced diameter to insure the proper location of the junction. For 
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differential thermocouple readings, the junction must also be elec- 
trically insulated from the blade unless the blade is used as one 
of the elements of the thermocouple. The insulation should be thin 
to reduce the thermocouple time lag and it should be of such a nature 
that the location of the junction is also known. In general, errors 
resulting from improper thermocouple location and installation can 
be reduced by placing the thermocouples as far apart as possible 
consistent with other requirements of the installation and by placing 
the thermocouples in small shallow holes. • 


Effect of Thermocouple Hole on Heat-Flow Path 

An analysis has not been made to evaluate the effects of holes 
placed in the heat path similar to those holes in which thermocouples 
are placed. If the holes are small compared to the wall thickness 
and if the metal beyond the thermocouple junction is solid, the 
errors caused by the holes are believed to be negligible. 


Accuracy of Temperature Measurement 

Accurate measurements of the temperature difference between 
two known locations in the blade metal are necessary. In most cases 
the difference between the two observed absolute temperature readings 
is not accurate enough. The method of measurement that has been 
found to be most successful is the use of a differential thermocouple 
, circuit and the measurement of the potential from the circuit with a 
sensitive potentiometer and a light -beam galvanometer. Frequently, 
use of the blade metal as one of the thermocouple elements is con- 
venient because differential readings can be made with the thermo- 
couple wire bonded to the blade metal. The single-wire thermocouples 
can also be placed in smaller holes. In order to use this method of 
temperature-difference measurement, an accurate calibration of tem- 
perature against electromotive force must be made of the metals used 
to form the thermocouple. 

The necessity for the high degree of accuracy required in the 
differential temperature measurements between two locations in the 
blade metal can be illustrated by equation (la): 
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The outside heat- transfer coefficient is almost a linear func- 
tion of r for the range encountered in most investigations so 
that the coefficient is directly proportional to t x - ty and 
inversely proportional to t e - t x . The temperature difference 
t x - ty can range from less than 1° to over 30° F depending on the 
blade-metal thermal conductivity and the distance between the loca- 
tions where the measurements are made. A 1° error in reading this 
temperature difference can therefore amount to an error of from 
3 to 100 percent in the heat- transfer coefficient. The temperature 
difference t e - t x may range from 100° to 1000° F or higher, how- 
ever, so that considerably larger errors can be tolerated in its 
measurement. 


EXPERIMENTAL RESULTS 

Preliminary heat-transfer runs were made to determine how suc- 
cessful these methods of analyzing heat-transfer data are in actual 
practice. For the sake of simplicity, a single symmetrical water- 
cooled blade having a chord and span of 6 inches was mounted in a 
test section with contoured walls as shown in figure 4. With this 
wall arrangement, the pressure distribution found on a typical reac- 
tion turbine blade can be simulated. A series of heat-transfer runs 
were made with heated air at temperatures ranging from 200° to 
1000° F. 

The leading edge, the trailing edge, and the central sections 
of the test blade were of shapes that could be used with the 
temperature-distribution equations presented herein. Multiple ther- 
mocouples were placed in the line of heat flow at the leading and 
trailing edges and single thermocouples were placed in the wall at 
the central portion of the blade. The circles on the blade sketch 
indicate the locations of the temperature measurements. The blade 
used for these heat-transfer determinations was made of aluminum, 
which has a relatively high thermal conductivity; because of its 
high thermal conductivity, the data obtained from the single ther- 
mocouples in the blade walls could not be used to evaluate local 
heat-transfer coefficients at the central portion of the blade. 
Experimental data are therefore given for only the leading and 
trailing edges. 


Trailing-Edge Section 

Heat-transfer data obtained from the temperature measurements 
in the trailing edge are shown in figure 5. The scatter in the 
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experimental data points is quite small and the data fall in small 
clusters. Each cluster represents a different gas temperature. In 
this case, data were taken at 200°, 400°, 600°, 800°, and 1000° F. 
Variable coolant-flow points are shown at each gas temperature. The 
equation of the experimental line was determined by the method of 
least squares and the intercept Cl on the abscissa was 0.954. With 
this value of 0, the local effective gas temperature can be cal- 
culated at any of the observed gas temperatures in the series of 
runs represented on this plot. 

The slope of the line representing the data is ~r, which in 
this case is equal to -0.0477. A straight-line approximation of 
equation (25) for the trailing edge of this blade is 

h t Z 0.3775 T (25a) 

From this approximate relation, the resulting local value of the 
outside heat-transfer coefficient at the trailing -edge section is 
0.01802 (Btu/(°F)(sq ft)(sec)). The average outside coefficient 
for the entire blade was slightly higher than this value, which 
serves as a rough check on the magnitude of the trailing-edge 
coefficient. ' 

If the effective gas temperature is calculated from the local 
blade recovery factor and the local Mach number, the calculated heat- 
transfer coefficients are lower than for the case where the effec- 
tive gas temperature is calculated from the intercept on figure 5 . 

The calculation for the quantity of heat transferred is the same in 
either case, however, as long as corresponding values of heat- transfer 
coefficient and effective gas temperature are used in the calculation. 
Because it is usually easier to calculate the effective gas tempera- 
ture from a blade recovery factor, outside heat-transfer coefficients 
were also calculated at each gas temperature using this value of 
effective gas temperature. 

The recovery factor for this blade at the trailing edge was 
0.90 at a Mach number of 1.0. The Mach number was maintained con- 
stant at 1.0 for this series of runs. Use of these values in equa- 
tion (52a) results in 


t e = 0.983 T (52b) 

A plot of outside heat-transfer coefficient against gas tem- 
perature where the effective gas temperature is defined by equa- 
tion (52b) is shown in figure 6. The data points shown were obtained 
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from arithmetically averaged blade temperatures at each gas temper- 
ature. All the experimental data used in figure 5 were also used 
to obtain figure 6. The outside heat-transfer coefficient increases 
as the gas temperature increases and is from 2.5 to 16 percent less 
than the constant coefficient obtained from the plot on figure 5. 


Leading-Edge Section 

Heat- transfer data obtained from the temperature measurements 
in the leading-edge section of the blade are shown in figure 7 for 
gas temperatures from 200° to 800° E. The scatter in the data for 
the leading-edge section of the blade is considerably greater than 
that for the trailing-edge section (fig. 5). There are two explana- 
tions for this scatter: (l) The change in local heat- transfer coef- 

ficient with distance from the stagnation point is very rapid. In 
this blade, the heat-transfer coefficient was not measured directly 
at the stagnation point and any stream fluctuations that might cause 
the stagnation point to shift would materially affect the value of 
the local outside heat-transf er coefficient. (2) In addition to 
this condition, the coolant flow near the leading edge of the blade 
was unstable and caused the blade temperatures to be unstable also. 
In a previous discussion, small errors or fluctuations in the blade- 
temperature measurements were shown to affect the accuracy of the 
data; the effect is illustrated here. 

At the leading edge of the blade, the value of £2 was found 
to be 0.940. The straight-line approximation of equation (20) for 
the leading-edge section of this blade is 

h t ~ 1.738 T (20a) 

By using this relation, the local heat-transfer coefficient at the 
leading-edge section was found to be 0.0441 (Btu/(°F)(sq ft) (sec)). 

Theory for the leading edge of cylinders indicates that the 
heat-transfer coefficient at the leading edge should increase with 
increase in gas temperature for a constant weight rate of gas flow. 
The plot of heat-transfer data in figure 7 indicates that this trend 
is also present here because the data points could be best repre- 
sented by a curved line. The data obtained at a gas temperature of 
1000° F were not included in figure 7 because of this trend. The 
data points at gas temperatures from 200° to 800° F can be approxi- 
mately represented by a straight line and were presented in this 
figure for illustrative purposes only. 
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For a more accurate determination of the heat-transfer coeffi- 
cients, the leading-edge- section data were also calculated using 
effective gas temperatures calculated from the recovery factor and 
the Mach number. The Mach number at the leading edge was 0.37 and 
the local recovery factor was 0.87, so that from equation (52a) 

t e = 0.997 T (52c)' 

A plot of heat-transfer coefficient against gas temperature at 
the leading edge of the blade where the effective gas temperature 
is defined by equation (52c) is shown in figure 8; for reference 
the heat- transfer coefficient obtained from figure 7 is shown as 
a dashed line. 

The heat- transfer coefficient at the leading edge is dependent 
to a very large degree on the leading- edge configuration. The coef- 
ficient is lowest for a circular section, increases materially as 
the section becomes elliptical, and becomes infinite for the leading 
edge of a flat plate. The heat-transfer coefficients represented 
by the solid line on figure 8 lie between those calculated from 
theory for a circular section and an elliptical section with a major- 
to minor-axis ratio of 2. Although this fact cannot be used as an 
accurate check on the heat-transfer data at the leading edge, it 
shows that the data are in the right range and the trend of varia- 
tion in gas temperature is verified. 


CONCLUDING REMARKS 

The analytical methods for the determination of local heat- 
transfer coefficients and effective gas temperatures present a con- 
venient method for obtaining these data and the actual application 
of the equations to experimental data is relatively simple. The 
application of the equations for the various sections of the blade 
is made in a similar manner for each section, greatly increasing 
their utility. 

Measuring the temperature gradient in the metal with two or 
more thermocouples whenever possible is advantageous because the 
analysis methods using one thermocouple are considered less accurate 
and considerably more difficult. 
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A desirable feature of these analytical methods is that they 
are not limited to turbine blades alone but can be used to determine 
local heat-transfer coefficients and effective temperatures for many 
other types of apparatus where the heat is being transferred through 
metal in a shape that can be approximated by the shapes discussed 
herein. 


Levis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 10, 1950. 
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APPENDIX A 
SYMBOLS 


A 

B 

C 

C P 

d 

E 

F 

& 

f 

G 

% 

E 1 

h 

I 

i 

J 0 

Jl 


The following symbols are used in this report: 
blade surface area, sq ft 



kg sin a 
constant 

specific heat at constant pressure, Btu/(lb)(°F) 
differential 

quantity in equation for trailing-edge analysis 
(equation (30)) 

radiation geometry factor 

radiation geometry factor for gray surfaces 
function 

quantity in equation for trailing-edge analysis 
(equation (46)) 

Hankel function of zero order 

Hankel function of first order 

heat- transfer coefficient, Btu/(°F)(sq ft)(sec) 

quantity in equation for trailing-edge analysis 
(equation (29)) 

'T'i 

Bessel function of zero order 
Bessel function of first order 
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k 

L 

M 

N 


n 

P 

P 

Q 

q 

r 

S 

T 
t 
w 
' X' 


y 

y 

Z,Z' 

a 

r 


constants of integration (function of ht/kg) 
thermal conductivity, Btu/(°F) (ft) (sec ) 
length, ft 
Mach number 

quantity in equation for trailing-edge analysis 
(equation (28)) 

exponent 

total pressure, lb/sq ft absolute 
static pressure, lb/sq ft absolute 
heat-flow rate, Btu/sec 

heat-flow rate per unit length, Btu/(ft) (sec) 
radius, ft 

quantity in equation for trailing-edge analysis 
(equation (46)) 

total temperature, °R 

temperature , °R 

weight-flow rate, lb/sec 

constant of integration (function of ht/kg) 
constant of integration (function of ht/kg) 
distance from trailing edge, ft (figs. 1(d) to (f)) 

i 

constants of integration (functions of h^/kg) 
trailing-edge wedge angle , deg 

function of ht/kg, thermocouple location, and blade 
configuration 
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r' 

7 

6 

6 

e 

t 

*1 

e 

A 

e 

^b 

T 

cp 

0 


intercept on plot of (1- T)/v c n against T 
(fig. 3(b)) 

ratio of specific heats 

thermocouple location dimension, ft (fig. 1(a)) 
blade-wall thickness, ft (fig. 1(a)) 
emiss'ivity 

fj_ (h t /kg, k B , and blade dimensions) 
fg (h t /kg, kg, and blade dimensions) 

angle, radians 
recovery coefficient 

thickness, ft (figs. 1(d) to (f)) 



Subscripts : 

av average 

B blade 

trapezoidal portion of trailing-edge section composed 
of trapezoid and rectangle (fig. 1(e)) 


b 
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b' 


c 

d 


d' 


e 


g 

i 

o 

r 

T 

t 

w 

x 

y 


1 , 2 , 


first trapezoidal portion of trailing-edge section 
composed of tvo trapezoids (fig. 1(f)) 

coolant 

rectangular portion of trailing-edge section composed 
of trapezoid and rectangle (fig. 1(e)) 

second trapezoidal portion of trailing-edge section 
composed of two trapezoids (fig. 1(f)) 

effective; subscript used with symbol for temperature 
to denote temperature effecting heat transfer 

gas 

inside blade surface 
outside blade surface 
radiation 
thermocouple 

combination of radiation and convection 

radiating surface 

local value at some location 

local value at some location different from location 
of x 


7 refer to thermocouple locations unless otherwise noted 
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APPENDIX B 

DERIVATION OF HEAT-TRANSFER EQUATIONS 
Central Section of Blade 

The heat flow through the blade wall is assumed normal to the 
blade surface and the surface curvature small enough that the wall 
can be treated as a flat plate. By reference to figure 1(b) for 
symbols and dimensions, the following heat balance can be written: 

^ = h-t(t e - t 0 ) = (t Q - t±) = g- (tg - ti) = hi(ti - t c ) 

(Bl) 

where 

t Q metal temperature on outside surface of blade 
t^ metal temperature on inside surface of blade 

t B metal temperature at distance 6 from coolant passage 

The following equations can be obtained from equation (Bl): 

(B2) 
(B3) 


k B t i + ^e 
k B + 6h t 

k B t 0 -fc 5h i t c 
kg + 6hi 
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If r is defined here as 

p _ h t (k B + ShjJ 

k B^t + h i) + &h t h i 


then 

t B - t c = r(t e - t c ) 


( 12 ) 


(5) 


The effective gas temperature t e in equation (5 ) is defined 
as the temperature that the test blade would assume if there were 
no heat transfer to or from the blade (reference 3). A blade recov- 
ery factor A b that is a function of the blade configuration and 
Mach number can be obtained, so that the effective gas temperature 
can be determined from the equation 


t e = T 


l + Ab 


7 - 1 


1 + 


7 - 1 



(B6) 


Similarly, the gas-temperature reading from a shielded ther- 
mocouple placed in the gas stream ahead of the blade will be 



(B7) 


where 


At 

■*r 


observed gas temperature reading from thermocouple, °R 

o 

total gas temperature, E 
recovery factor for thermocouple 
Mach number at thermocouple 


If 


*e 



(4a) 
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where Q is a constant of proportionality, from equations (B6), 
(B7), and (4a), 



In order to determine the effect of variations in local Mach 
numbers and recovery factors of the blade and the thermocouple, 
values of Cl were calculated for Mach numbers from 0 to 1.0 by 
use of values of recovery factors from 0.65 to 1.0. The maximum 
variation in 0 with Mach number was about 2-i percent. 

Combination of equations (5) and (4a) and division by tg yield 



For the case when two temperature measurements can be made to 
obtain the temperature gradient in the blade wall, as shown on 
figure 1(a), temperatures t]_ and tg at corresponding distances 

and 6 2 can be substituted into equation (B5), which results 
in two equations. Subtraction of the equation involving tg from 
the equation involving t-^ results in 


t l ~ t 2 _ h t h j( 5 l ~ 6 2 ) 

t e ~ t c k B (h t + hj_) + 'Sh t h i 

Equation (B6) can also be written 


(B9) 


t l " t c _ h t (k B + 5^) 

^e ” ^c kg(h^. + h^) + &h^.h i 


(BIO) 


Subtraction of both sides of equation (BIO) from unity and com- 
bination with equation (B9) to eliminate h^ give 


ti - t 2 



t l^ 


(Bll) 
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If r is now defined as 


^ < 5 1 ‘ 6 2 > 
r = — 

1 + ir (6 " & l) 
k b 

then equation (Bll) becomes 

tj. - t 2 = r(t e - t x ) 

Equation (B5) can be written 


(9) 


(7) 


where 


and 



(2) 


(16) 


(17) 


If t B in equation (2) is replaced by tg, equation (16) 
becomes 



(10) 


Leading-Edge Section of Blade 

Heat is assumed, transmitted from the gas stream to the blade 
coolant at the leading edge of the blade along a sector with an 
included angle d 9, as shown in the leading-edge cross-sectional 
view in figure 1(c). It is assumed that no heat is transmitted in 
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a direction normal to the sector path from the gas to the coolant; . 
this assumption has been verified by relaxation calculations for 
well-rounded leading-edge sections. The following heat-balance 
equations can now be written: 


= (*e - t o)( r o 

d0)h t 

(B12) 

dq = (t ± - t c )(r 1 

d0)hi 

(B13) 

dt B 

dq = r d0 k B — 


(B14) 


where 

dq heat transferred along sector having included angle d0 

per unit length normal to sector, Btu/(ft) (sec) 

r, r.^, r Q radii to concentric arcs of circles 

t B blade temperature at any point in sector, °E 

t i inside-blade-surface temperature, °E 

t Q outside-blade- surface temperature, °E 

From equation (B14), 

t B-kffe log « C r < B15 > 

Then 

*o Cr 0 (B16) 

and 

t i-^B loe * Cr i <B17) 

Subtracting equation (B17) from (B16) yields 

h-p d0 , . 

dq = S (to - ti) (B18) 

. r o 
lo Se — 

1 i 
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It is obvious that 


4 o - *i 


t D - H 


t e - t c (t e - t Q ) + (tj_ - t c ) + (t Q - t ±l 


(B19) 


Combination of equations (B12), (B13), (B18), and (B19) yields 


t - t. = 
o 1 


te " t c 

k B / 1 ( 1 

, r o W r o + h i r i ( 
logp — N 


(B20) 


Subtraction of equation (B17) from (B15) and division by (B18) give 


% ~ t i 
t o ~ t i 


lo *e^ 


(E21) 


When equations (B13) and (B18) are combined, 


t i - *c = 


r i h i lo ®e ~ 


— (*o " ti> 


(B22) 


Combination of equations (B20) and (B22) results in 


ti - t c - 


k B te ~ tc 

, , _ r o , k B ( 1 1 

r-h. log — 1 + — It— — + t— - 

11 r i log v^o k i r i 


(B23) 


From a combination of equations (B20) and (B2l), it is found that 


loge r- 


lo Se — 1 + 


t e - t c 


lo Se t : 


(B24) 
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It is obvious that 


t B " t c = (*B ‘ ti) + (*i " t c) 
By combination of equations (B23) to (B25), 


(B25) 


■'B 


log„ — + , 

&e r 4 h,- r 


B 


1 i 


~ t c 


log e -2 + 


V B 


r i 


h t r o 


h i r i 


(B26) 


If two temperatures tj and t 4 are measured in the sector 
at points 3 and 4 at radii r 3 and r 4 , these values can be sub- 
stituted into equation (B26), which results in two equations - one 
for each temperature. Subtracting the equation involving t 4 from 
the equation involving t-j yields 


s - u . 

- ‘c 


r h, r_ 

o t .. 3 

log e — 


‘■B 


r 4 


\ k B 

Equation (B26) can also be written 


1 + h t log e — + 


h i r i / 


(B27) 


*5 " t c 
t e - t c 


ht ^ loge ^ + ^ 
VB h l r i/ 

i , t. f r o i r o . r o 

1 + h tl 7— l°Se — + 


C B 


h i r i. 


(B28) 


Subtraction of both sides of equation (B28) from unity and combina- 
tion with equation (B27) give 


t 4 =1 


r o h t w r 3 

IT - 1 Se ~ 

*B 4 

. r o h t - r o 
d + — logg — 


“ ^3 ) 


(B29) 


K B 


r 3' 
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If T is defined as 


r = 


r h + r, 

o t . 3 

— log — 

k B ® r 4 


1 + 


r o h t 

k B 


( 20 ) 


log, 


e T? 


equation (B29) can be written 

t 3 - t 4 = T(t e - t 3 ) (18) 

Combination of equations (18) and (4a) and division by tg yield 


= r n . 


b g 




(19) 


In order to determine the value of the inside heat-transfer 
coefficient hj_, both sides of equation (B26) are subtracted from 
unity and the equation is rearranged to yield 


h-t = 


r i L r o ^ k B 
k B l log e r h t r 0j 


^B 


J B/ 


r i L r o + 

’ k B V° ge r i h t r o. 


If i and tj are defined as 


(B30) 


k B l h t r o 


t =ii(^ + iog e ^) 


(21) 


and 


1 

*1 = n 


B 


k B \ h t r o 
equation (B30) can be written 


+ log e — 


( 22 ) 


hi = 


£ (M>V* 


“B/ 


( 2 ) 
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Trailing-Edge Section of Blade 

The trailing- edge cross section of most turbine blades can be 
approximated by a trapezoid, a rectangle, or a combination of 
trapezoids and rectangles. Heat-transfer equations are given for 
each case. 


Trapezoidal trailing- edge section. - The temperature distribu- 
tion in a trapezoidal trailing-edge section of a turbine blade is 
derived in reference 4 and is given herein by the expression 



h.N 

i 

I + h t E 


(B31) 


vhere 


N = 


t2 


2B^k 


B 


J()(i!) + iEqCiI) 


I = iJ^(i|-^) %(i£g) “ %(i£l) 

E = - U(i^l) JoUtg) + 'iJi(ili) iHo(iS 2 )l 

2B^k 0 L - 1 


i = 2B 


it _ \ _ 7 1 - tan a 

^ + 4 1 / 1 \ 2 tan a . 


B 2 = 


kg sin a 


a = tan' 


•i T 5 " T 1 
2L 


ii 

^2 

y > > "^3 > ^ 


(28) 

(29) 

(30) 


evaluated for y + j = 0 
evaluated for y = L 

shown on trailing-edge section for trapezoidal shape 
(fig. 1(d)) 
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correction term in equation for £ to account for 

rounded end of wedge and considered more accurate than 
T^/2 as given in reference 4. 

If two temperatures, tg and tg, are measured in the trailing 
edge at distances yg and yg, as shown on figure 1 (d), these 

values can be substituted into equation (B31), which results in two 
equations, one for each temperature. Subtraction of the equation 
involving tg from the equation involving tg leaves 


- T 
4 T 1 





h i 


(B32) 


where N 5 and Ng are evaluated at the positions where tg and 
t 6 are measured. ' Eeplacement of tg by tg and N by Ng in 
equation (B31) and combination with equation (B32) result in 


t 5 



* %) 


If r is defined as 


(B33) 


r = 



equation (B33) becomes 

t 5 " t 6 = ‘ 


(25) 


(23 ) 


Division of both sides of equation (23) by tg and combination of 
the result with equation (4a) yield 


t 


5 


- t 6 


g 




t 


(24) 
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Solving equation (B3l) for yields 


1 W 7 ^e ~ ^ E 

1 W - W 1 


(B34) 


* = N 
^ I 


* = Y 


results in a transformation of equation (B34) to 


h i " ft e - t c \ (2) 

Rectangular trailing-edge section. - The temperature distribu- 
tion in a rectangular trailing-edge section is given in reference 4. 
In the notation of this report, the expression becomes 

!e ; V I cosh ^ l! T j) (B35) 

6 C <p sinh cp^Lg + ^ + cosh cp^I^ + 5 


shown on rectangular portion of trailing-edge section 
sketch (fig. 1(e)) 


correction for semicircular end of trailing-edge section 
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#> 


and 


cp = 



If two temperatures tg and tq are measured in the trailing 
edge at distances yg and y7, as shown on figure 1(e), these 

values can be substituted into equation (B35), resulting in two 
equations - one for each temperature. Subtraction of the equation 
involving tg from the equation involving tj leaves 



hi 

cosh Cp| 

v y 7 + 1 n) 

| - cosh cp(y 5 + | Tl )_ 

°Pk E 

, sinh Cp| 

( L 2 + f T l, 

) + ^ cosh cp(L 2 + | T 1 ) 


(B36) 


Beplacement of tg by tg and y by yg in equation (B35) and 
combination with equation (B36) result in 


- t7 . (te _ t5)(B37) 

cosh tp(y 5 + 2 


By letting 


cosh cp^y ? + ^ - coshcp^y 5 + | 


cosh cp 


( y s + f T 0 


(33) 


equation (B37) becomes 


t 5 - t? - T(t e - t 5 ) 


(31) 


Division of equation (31) by tg and combination of the result with 
equation (4a) yield 


tc - t 7 / t q N 

-hr^ - r n - r 


( 32 ) 
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Equation (B35) can be solved for h^, which results in 


h i = 


kgCP 


By letting 


cosh cp(y + f T x ) ( t e - t \ 
> sinhCp(L 2 + J T^) Ve ~ t B/ 


cosh cp(L 2 4 | Tl ) 
k B cp sinh cp(^L 2 + f 

( B38 ) 


t = 


cosh cp 




and 


B = 


k# sinh < 5 p(lig + f T X ) 

cosh cp^L 2 + | T X ) 
kgCp sinh cp + “ T"' 


(34) 


(35) 


equation (B38) becomes 


h i 


*e - V 
V t e " t B j 


( 2 ) 


- Tl 


Combination trapezoidal and rectangular trailing edge, - By 
writing equation (B35) for tg and tq~, equating the two expres- 
sions, and solving them, the temperature distribution in the rec- 
tangular portion of the trailing edge can be written as 


t B = (^e " *7^ 


cosh 


cosh cp 


Hr + l \) 

( l 2 + f T i) 


(B39) 


where 

tg temperature at any point in rectangular portion 
t y temperature at point where y = Lg = 7-j 
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The temperature distribution for any trapezoidal shape is given 
in reference 8, and in the notation of this report becomes 

t e - t B = KJ 0 (i£) + ZiBo(il) (B40; 

The integration constants K and Z are functions of hj_ . 

At the junction of the two portions of the trailing edge, the metal 
temperature and the heat flow must be continuous; therefore relations 
between the integration constants can be written. 

At the junction, the blade temperature is t^; therefore 

t e " t? = KJ 0 (i£ b>1 ) + ZiH 0 (i| bA ) (B41) 

where 

lb i evaluated for y = 0 
and 

b = 

l T 3 - T 1 
% = tan -1 — - 

For the heat flow to be continuous, the following relation must exist: 

dt b dt d 

KBftb TJ = k B»d (M2) 

but A b = A d at the junction where the subscripts b and d refer 
to the portion of the trailing-edge section in figure 1(e). 

Differentiating equations (B39) and (B40) with respect to y 
and substituting the result in equation (B42) result in 

(t e - t7 ) cp£ tanhCp d ^L 2 + J Tl) = [- KiJi(i|b,l) + ZHi(i£b,i)j 



(B43) 
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Temperatures tg and at corresponding distances y B 

and y 7 can be substituted into equation (B39), which results in 
two equations. Subtracting the equation involving tg from the 
equation involving t -7 results in 

cosh cp d (y 7 - + | tA - cosh cp d (y 5 + f T X ) 

t 5 - t 7 * — - — * — L (t e - t 7 ) 

cosh | 

(B44) 


Equation (B39) can also be written 


t 


e 



cosh q> a (y 5 + f Tl ) 
cosh + 2 ^ j 


(B45) 


The rigorous solution to equations (B44) and (B45) is rather 
complex. From equations (23) and (31) it can be seen that the final 
solutions for the various shapes take on the general form 


t 


x 




( 1 ) 


where T is a function of the blade dimensions , the end conditions, 
and the ratio h^/kg. 


The simplest method of solving equations (B44) and (B45) to get 
them in the general form of equation (l) is to calculate the slope 
T d of a plot of tg - t 7 against tg. In order for the slope to 

be constant it cannot involve h d . The slope of the line is 


^(^5 " ^7) 






d(tc ) 


dhi 


(B46) 


From equation (B44), 


d(tc- - t 7 ) 


dh i 


cosh cp a (y 7 + f T x ) - cosh 4 - | T.J dt. 


cosh cp d (L 2 + -d T- 


n 

4 


dhi 


(B47) 
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and from equation (B45), 


dt 5 _ cosh ^( 7 5 + f T l) dt 7 
cosh cp^l^ + X T l) . dili 

Then from equations (B46) to (B48), 

r _ cosh <P a (y 7 + f T i) - cosh <p d (y 5 ♦ f Tj) 

COSh CP d (y 5 + J Tl ) 


(B48) 


(33a) 


For the trapezoidal portion of the trailing-edge section in 
figure 1 (e), the temperatures t^ and tg can he substituted into 

equation (B41), which results in two equations. Subtracting one 
equation from the other and differentiating yield 




dh ± 


M 

dh 


^o( i ^b,6) “ Jo^^b,?)] + iH o( i ^b,6) -iHoCitb,?) 

(B49) 


and 


dty 

dh[ 


lb JQ^^b,?) - -§b 


(B50) 


where 

^b ,6 evaluated for y = y 6 

^b ,7 evaluated for y = y-j (y 7 = 0 for this case) 

Division of equation (B49) by equation (B50) to obtain the value 
of the slope for the trapezoidal part of the trailing edge 

results in 


V 


J C>( i £b,6) ' J 0^^b,7) + if " iH 0^^b,7^ 


Jo^b,?) + || iHo(i£b,7) 


(38) 
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The value of dz/dK in equation (38) can be found from equa- 
tions (B4l) and (B43) as follows: 

When equations (B41) and (B43) are differentiated with respect 
to h^. 


dt 7 _ dK 
dh d - " dhj_ 




Si 


(B51) 


and 


dt 

dh 


\ ^d tanh <Pd( L 2 + f T l) = 


2 V 

b,l 


fr iJi(itb,i) 


Sr 

(B52) 


where 

Ib^i evaluated for y = 0 
Elimination of dt 7 /dh i yields 


2B b iJ l^ i ^b,l) + ^b,l %. 

tanh cp d 

( L 2 + f T l) 

J 0 ^ 1 ^b , 1 ) 

2B b M^b^l) ' £b,l ^ajjtanh c Pd( 

L 2 + f T l)_ 

iSo^^b,!) 


(39) 

Combination of two trapezoidal sections in trailing edge. - 
The temperature distribution in the trapezoidal portion b' is given 
in different form in reference 8 as 

t e - t B = K'J 0 (iS b .) + Z'i%(iS b ,) (B53) 

and in portion d f the distribution is 

t e - tg = X'Jo(itd') + Y'iHoUtdO (B54) 


where 

K'jZ’jX^Y' integration constants 
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and • 


i 


B d i = 1 


k B sin a^» 


Here again the integration constants are functions of h^ and 
at the junction of the two trapezoidal sections the metal temperature 
and heat flow must be continuous. The following relations between 
the integration constants can be written: 


For temperature continuity, 

' K ' J 0^V,1> + Z ’ iH 0 (1 fbM> - + Y ' 1H 0 (l5 d',2> 

(B55) 

and for heat-flow continuity. 


, . dt b’ , n dt d- 
k B A b' = k B A d' 


(B42a) 


Differentiating equations (B53) and (B54) with respect to y 
and substituting in equation (B42a) yield 
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r-^“ [- x ) + -l)] 

^ b ',1 L J 

2 

= T^~g ^ x ' iJ l( i ^d' , 2 ^ + Y ' H l( i ^d',2^ (B56) 

where 

A b* = A d' 

^b» 1 evaluated for y = 0 

^d' 2 evaluated for y = Lg 

The equation for portion d' is 

ts - t? = rd> (t e - tg ) (40) 

From equation (38), the equation for T^i can be written by 
inspection: 


r*. = 


J (il )-j (it ) + gl("iH(i£ ) - iH (i| ) 

0 cL',7 0 cl',5 dX' 0 d*,7 0 d',5_ 


J 0^d',5) + i ®o( i £d' ,5^ 


(B57) 


where 

^d ’ ,7 evaluated for y = y^ 
ld » ,5 evaluated for y = y 5 

Y',X' integration constants for portion d' that replace integra- 
tion constants K and Z for portion b in equation (38) 
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The equation for dY'/dX' can also be written by inspection 
from equation (39). If the rectangular portion d were removed 
from portion b (fig. 1(e)), the quantity (Lg + — T^) in equa- 
tion (39) would be equal to zero and the equation would reduce to 


dZ 

dX 


iJ i (l£ b,i> 



(39a) 


Because nothing is attached on the rear of portion d' on 
figure 1(f), the equation for dY’/dX’ is 


dY' 
dX' 

where 

. It 

, evaluated for y + — T = 0 
d' ,1 ^ 4 1 

and 


iJ l^d',l) 

Hl^d* ,l) 


(B58) 


^d' _ 2B d' 

A correction for the rounded end is included in equation (B58) 
but not in equation (39a). 

Substitution of equation (B58) into (B57)-and notation of the 
definition of N given in equation (28) yield 




(42) 


When the same process as before is used to obtain an equation 
for portion b' in the form 


t 7 _ t 6 = (t e - t 7 ) 


(43) 
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it is found that 


rv, = 


- ~ J Q ^ b 1 , 7 ^ ' ^O^by/j] 


J o( i ^h',7) + 7§7 ^ 1 , 7 ^ 


The value of cLZ'/dK* can be obtained from equations (B55) and 
(B56) by differentiating and grouping to yield 

J 0^^b',l) + ^7 iH o( i tb , ,l )= ^l|5o (i ^d , ,2) + §T iBo(iSd* , 2 >] 


(B59) 


IT iJ l^b',l) + ft? H l(^b',l) 


s[- 


dY’ „ n B d'^bM 


(^d' ,2) + d^T 


n s d'^M 

> 2 _L 2 t 

- JB b' ^d',1 


(B60) 


Elimination of dX'/dh^ from equations (B59) and (B60) and 
substitution of the value of dY'/dX' from equation (B58) yield 

dZ* = G1J 1 ^^m) + SJ 0^ 1 ^b',l^ / 46 

dK 1 " GH^iev^) - SiHo(i| b . #1 ) 


G = J 0 (Ud* ,2) + 


1J l(i^d l ,l) ^(^d* ,2) 


B d ,2 ^b' ,1 r iJ l^d‘ ,l) H 1 C d* , 2 ) .w. t > 

TW~^ HiUU .) ’ lj i (l5 a', 2 ) 

B b' M',1L l v s d*,l' J 
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(b) Plot II. 

Figure 3. - Illustrative graphical solution of 
temperature-distribution equation for central 
portion of blade using one thermocouple. 
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Figure 4. - Water-cooled "blade used for experimental determination of local heat-transfer 

coefficients and effective gas temperatures. 















= 0.997 T 
, = 0.940 T 
rom fig. 7 ) 
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•e 8. - Comparison of outside-surface heat- transfer coefficients at leading edge 
blade obtained from figure 7 with those calculated where effective gas tempera- 
•e was based on recovery factor. 



